Enhancing heat transfer in nanofluids by carbon nanofins 
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Nanofluids are suspensions of nanoparticles and fibers which have recently attracted much atten- 
tion due to their superior thermal properties. Here, nanofluids are studied in the sense of nanofins 
transversally attached to a surface, so that dispersion within a fluid is mainly dictated by design and 
manufacturing processes. We focus on single carbon nanotubes thought as nanofins to enhance heat 
transfer between a surface and a fluid in contact with it. To this end, we first investigate the ther- 
mal conductivity of those nanostructures by means of classical non-equilibrium molecular dynamics 
simulations. Next, thermal conductance at the interface between a single wall carbon nanotube 
(nanofin) and water molecules is assessed by means of both steady-state and transient numerical 
| experiments. Numerical evidences suggest a pretty favorable thermal boundary conductance (order 

. of 10~ 7 [Wm~ 2 K~ x ]) which makes carbon nanotubes ideal candidates for constructing nanofinned 

^—i • surfaces. 

(N 

PACS numbers: 66.70.+f, 05.10.-a, 62.25.+g, 68.35.Md 

. I. BACKGROUND 

Nanofluids are suspensions of nanometer-sized solid particles and fibers, which have recently become a subject 
of growing scientific interest because of reports of greatly enhanced thermal properties P, 0] ■ Filler dispersed in a 
nanofiuid is typically of nanometer size, and it has been shown that such nanoparticles dispersed in a base fluid are 
able to endow it with a much higher effective thermal conductivity than pure fluid 0, H| : Significantly higher than 
, those of commercial coolants such as water and ethylene glycol. In addition, nanofluids show an enhanced thermal 
conductivity compared to theoretical predictions based on the Maxwell equation for a well-dispersed particulate 
composite model. These features are highly desirable for applications, and nanofluids may be a strong candidate 
for new generation of coolants 0. The use of nanofluids in many industrial sectors, including energy supply and 
production, transportation and electronics appears promising. A review about experimental and theoretical results 
£Sj , on the mechanism of heat transfer in nanofluids can be found in Ref. [f| , where Authors discuss issues related to the 
' technology of nanofiuid production, experimental equipment, and features of measurement methods. A large degree 
of randomness and scatter has been observed in the experimental data published in the open literature. Given the 
inconsistency in these data, it is impossible to develop a comprehensive physical-based model that can predict all the 
trends. This also points out the need for a systematic approach in both experimental and theoretical studies @. 

In particular, carbon nanotubes (CNTs) have attracted great interest for nanofiuid applications, because of the 
claims about their exceptionally high thermal conductivity [7(. However, recent experimental findings about CNTs 
report an anomalously wide range of enhancement values that continue to perplex the research community and 
remain unexplained [8|. For example, some experimental studies showed that there is a modest improvement in 
thermal conductivity of water at a high loading of multi- walled carbon nanotubes (MW-CNTs), ~35% increase for 
a 1 wt% MWNT nanofiuid [|. These authors attribute the increase to the formation of a nanotube network with 
a higher thermal conductivity. On the contrary, at low nanotube content, <0.03 wt%, they observed a decrease in 
thermal conductivity upon an increase of nanotube concentration. On the other hand, more recent experimental 
investigations showed that the enhancement of thermal conductivity as compared to water is varying linearly when 
MW-CNT weight content is increasing from 0.01 to 3 wt%. For a MWNT weight content of 3 wt% the enhancement 
of thermal conductivity reaches 64% of the base fluid (e.g. water). The average length of the nanotubes appears to 
be a very sensitive parameter. The enhancement of thermal conductivity compared to water alone is enhanced when 
nanotube average length is increasing in the 0.5-5 /im range fipj . 



Clearly, there are difficulties in the experimental measurements [ll| . but the previous results also reveal some 



underlaying technological problems. First of all, the CNTs show some bundling or the formation of aggregates 
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originating from the fabrication step. Moreover it seems reasonable that CNTs encounter poor dispersibility and 
suspension durability due to the aggregation and surface hydrophobicity of CNTs as a nanofiuid filler. Therefore, the 
surface modification of CNTs or additional chemicals (surfactants) have been required for stable suspensions of CNTs, 
because the base fluid for the coolant has polar characteristics. In the case of surface modification of CNTs, water- 
dispersible CNTs have been extensively investigated forpotential applications, such as biological uses, nanodevices, 
novel precursors for chemical reagents, and nanofiuids Q. A popular solution for increasing dispersion of CNTs is 
based on functionalization. Oxygen-containing functional groups have been introduced on the CNT surfaces and 
more hydrophilic surfaces have been formed during this treatment, which enabled to make stable and homogeneous 
CNT nanofiuids [12]. Alternative solutions relay on ultrasonic disrupting, which significantly decreases the size of 
agglomerated particles and number of primary particles in a particle cluster, such that thermal conductivity increases 
with the elapsed ultrasonication time [131 ]. In any case, it is clear that many parameters affect the thermal conductivity 
including size, shape and source of nanotubes, surfactants, power of ultrasonic, time of ultrasonication, elapsed time 
after ultrasonication, pH, temperature, particle concentration and surfactant concentration [lij ]. Hence, there is a lot 
of room for technological optimization. 

In the above brief review, only the free suspensions of CNTs have been considered, i.e. nanofiuids were the 
highly thermally conductive filler is free to move. In fact, beyond the favorable aspect ratio, the CNTs dispersed 
in nanofiuids lead to enhanced thermal conductivity, which cannot be explained by traditional conductivity theories 
such as the Maxwell's mixing theory or other macroscale approaches. Recently, Jang & Choi [l5[ have found that 
the Brownian motion of nanoparticles at the molecular and nanoscale level is a key mechanism governing the thermal 
behavior of nanoparticle-fluid suspensions. They have devised a theoretical model that accounts for the fundamental 
role of dynamic nanoparticles in nanofiuids. Essentially, these authors discovered a fundamental difference between 
solid/solid composites and solid/liquid suspensions in size-dependent conductivity. In the original reference [l5j . they 
claim that, even though the random motion of nanoparticles is characterized by zero time average, the vigorous and 
relentless interactions between liquid molecules and nanoparticles at the molecular and nanoscale level translate into 
conduction at the macroscopic level, because there is no bulk flow. 

However there is another way to exploit the concept of nanofiuids, i.e. to modify the base fluid properties by 
nanostructures. Essentially one promising way consists in fixing the relative position of the CNTs on nanoengincercd 
surfaces. This means moving from nanotubes (where the focus is on the shape of the nanostructure) to nanofins 
(where the focus is on the function of the nanostructure). In this way, the Brownian motion of nanoparticles is 
gone, but there is still the possibility to enhance the heat transfer by increasing the effective cooling surface area. 
The fundamental technological advantage is that the dispersion of CNTs is controlled by design and only limited 
by the manufacturing process. Nowa day s efficient cooling of silicon chips using microfin structures made of aligned 
MW-CNT arrays has been achieved [la] . The tiny cooling elements mounted on the back side of the chips enable 
power dissipation from the heated chips on the level of modern electronic devices demands. Nanotubes utilized as 
thermal fins (nanofins) are mechanically superior compared to other materials being ten times lighter, flexible, and 
stiff at the same time [H| . Nanofins are extensively investigated also from a modeling standpoint [ijj • The current 
challenge is to develop industrial manufacturing processes for macroscopic growth of carbon nanotube mats fl8| . 

The present paper aims to investigate by molecular mechanics based on force fields (MMFF) the thermal perfor- 
mance of nanofins made of single wall CNTs (SW-CNTs) cooled by water. This work focuses on the astonishing 
thermal properties of these nanostructures, in particular, when they interact with the surrounding base fluid. The 
single wall CNTs were selected mainly on the basis of time constraints due to our parallel computational facilities. The 
following analysis can be split into two parts. First of all, the heat conductivity of SW-CNTs is estimated numerically 
by both simplified model (Section 1, where this approach is proved to be inadequate) and detailed three-dimensional 
model (Section 2). This first step is required for tuning the numerical model and validate the vacuum results with 
literature results. Next, the thermal boundary conductance between SW-CNT and water is computed by two meth- 
ods: the steady state method (Section 3.1), mimicking ideal cooling due to strong Brownian motion, and transient 
method (Section 3.2), taking into account only weak Brownian motion. 



II. HEAT CONDUCTIVITY OF SINGLE- WALL CARBON NANOTUBES: A SIMPLIFIED MODEL 

In order to significantly downgrade the difficulty of studying energy transport processes within a carbon nanotube, 
authors often resort to simplified low dimensional systems such as one-dimensional lattices [l9l - l23 |. In particular, heat 
transfer in a lattice is typically modeled by the vibrations of lattice particles interacting with the nearest neighbors 
and by a coupling with thermostats at different temperatures. The latter are the popular numerical experiments based 
on non-equilibrium molecular dynamics (NEMD). In this respect, to the end of measuring the thermal conductivity 
of (5, 5) single wall nanotubes (SWNT), we set up a model for solving the equations of motion of the particle chain 
pictorially reported in Fig. 1 where each particle represents a ring with 10 atoms in the real nanotube. In the present 
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model, Carbon-Carbon bonded interactions between first neigbors (i.e. atoms of the ith particles and atoms of the 
particles i ± 1) separated by a distance r are taken into account by a Morse- type potential expressed in terms of 
deviations x — r — rg from the bond length tq: 



V b (x) = V (e 



2e~ 



where Vq is the bond energy while a is assumed a = Tq/2. Following |2f 
distance between two consecutive particles at equilibrium is assumed ro = 
the total force, Fi, acting on the ith particle is computed as: 



(1) 



, bond energy is Vb = 4.93[eV], while the 
0.123[nm]. At any arbitrary configuration 



Ft 



N bon sin <d 



dV b dV b 

— {dxi-i) + — {dx i+1 ) 



(2) 



with dxi-j = x\ — Xi-j, dxi+j = Xi+j — Xi and N bon denoting the number of Carbon-Carbon bonds between two 
particles, whereas a penalization factor sin?9 can be included to account for bonds not aligned with the tube axis. 
Here, we use free-end boundary condition, hence forces experienced by particles at the ends of the chain read: 



N bon sin •& 



dv b N 

^ {dx2 ' 



F, 



N 



-N bon sin$ 



dV b 



(3) 



Let Pi and m; be the momentum and mass of the ith particle, respectively, the equations of motion for inner particles 
take the form: 



dxj _ Pj_ djh _ p 
dt m, ' dt 



(4) 



whereas the outermost particles (i = 1, N) are coupled to Nose- Hoover thermostats and are governed by the equations: 



dxj _ Pi dpi _ 
dt m , ' dt 



KPl ' dt Q 



2m, 



= ^ 



(5) 



with fcf,, To, Nf and tt denoting the Boltzmann constant, the thermostat temperature, number of degrees of freedom 
and relaxation time, respectively, while the auxiliary variable £ is typically referred to as friction coefficient 1261 . 
Nose-Hoover thermostatting is preferred since it is deterministic and it preserves canonical ensemble (see, e.g., [27] 
and [28j for further details on thermostats in molecular dynamics simulations). 

Local temperature Tj(i) at a time instant t is computed for each particle i using energy equipartition: 



Ti(t) 



1 



k b N f 



(6) 



where () denotes time averaging. On the other hand, local heat flux Ja transferred between particle % and i + can 
be linked to mechanical quantities by the following relationship [U [23| : 



Ji = ( tt- (dxi+i) 

, mi ox 



(7) 



The above simplified model has been tested in a range of low temperature (300[!f] < T < 1000[K]), where we 
noticed that it is not suitable to predict normal heat conduction (Fourier's law). In other words, at steady state 
(i.e. when heat flux is uniform along the chain and constant in time), it is observed a finite heat flux although no 
meaningful temperature gradient could be established along the chain (see Fig. 1). Thus, the above results predict 
a divergent heat conductivity. Here, it is worth stressing that one-dimensional lattices with harmonic potentials are 
known to violate Fourier's law and exhibit a flat temperature profile and divergent heat conductivity. On the other 
hand, consistently with the present numerical experiments, it has been demonstrated that anharmonicity alone is 
insufficient to ensure normal heat conduction 11911 - 



III. HEAT CONDUCTIVITY OF SINGLE-WALL CARBON NANOTUBES: DETAILED THREE 

DIMENSIONAL MODELS 

In all simulations below, we have adopted the open-source molecular dynamics (MD) simulation package GROningen 
MAchine for Chemical Simulations (GROMACS) [29M3l| in order to investigate the energy transport phenomena in 
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FIG. 1: One-dimensional model: Lattice chain of particles interacting according to a Morse- type potential {TJ. End-particles are 
coupled to Nose-Hoover thermostats at different temperature (Thot = 320[-K] and T co id = 280 [if]). Despite of the anharmonicity 
of the potential, normal heat conduction (Fourier's law) could not be established. Here, heat flux is computed by eq. ©. 
However, consistent results are obtained by eq. ()12p below which predicts: {(hot) k^Thot = — {(cold) kbT co id = 1.11 x 10 _7 [W]. 



three-dimensional SWNT obtained by a freely available structure generator (Tubcgen) [32j . Three harmonic terms 
are used to describe the Carbon-Carbon bonded interactions within the SWNT. Namely, a bond stretching potential 
(between two covalently bonded carbon atoms i and j at a distance r$j ■): 

V b (r ij ) = lkl j {r ij -r° j )\ (8) 

a bending angle potential (between the two pairs of covalently bonded carbon atoms and (j, k)) 

V a (6 ijk ) = i k% h (cos 6 ijk - cos 6% k ) 2 , (9) 

and the Rychacrt-Bcllcmans potential for proper dihedral angles (for carbon atoms i, j, k and I) 

1 

2 



V rb (<Pi 3 kl) = ^kf jM (1 - COs20y fei ) (10) 



are considered in the following MD simulations. Here, 9^ and <pijki represent all the possible bending and torsion 
angles, respectively, while r+j = 0.142[nm] and #°- fe = 120° are reference geometry parameters for graphene. Non- 
bonded Van der Waal interaction between two individual atoms i and j at a distance r,j can be also included in the 
model by a Lennard- Jones potential: 



Vnb = 4ecc 



/ \ 12 / \ 6 ' 
/ &C C \ ( vcc \ 



(11) 



where the force constants k\^ kf ik and kfj kl in ijHJ), ©, (fTTJ|) and parameters (dec, ecc) in (JXTJ) are chosen according to 
the table 1 below (see also [33| and [34|). In reversible processes, differentials of heat dQ rev are linked to differentials 
of a state function, entropy, ds through temperature: dQ rev = Tds. Moreover, following Hoover (2^,(3^1, entropy 
production of a Nose-Hoover thermostat is proportional to the time average of the friction coefficient (£) trough the 
Boltzmann constant kb hence, once a steady state temperature profile is established along the nanotube, the heat flux 
per unit area within the SWNT can be computed as: 

« = ~(0^f, (12) 

where the cross section Sa is defined as Sa = 27rr&, with b = 0.34[nm] denoting the Van der Waals thickness (see also 
[36[). Here, the use of formula (|12[) is particularly convenient since the quantity (£) can be readily extracted from the 
output files in GROMACS. 
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FIG. 2: Three-dimensional model: Nose-Hoover thermostats are coupled to the end atoms of a (5, 5) SWNT. (Top) Both bonded 
© © (|10[) and nonbonded interactions are considered. In a three-dimensional structure, harmonic bonded potentials do 
give rise to normal heat conduction. Temperature profiles for two lengths (5.5[nm] and 10 [nm]) are reported. (Down) Several 
setups have been tested where some of the interaction potentials (j8j), (j9j), (|10|l and (|11[) are dropped out. BADLJ: Vb, Van, Vrb 
and V n b are consedered. BAD: Vb, V an , V r b are considered. BA: Vb and Van are considered. Bw denotes that Vb is computed 
with a smaller force constant k\j = 42000[fc Jmol~ 1 nm~ 2 ] according to [25j . 



The measure of both the slope of temperature profile along the inner rings of SWNT in Fig. 2 and Fig. 3 and heat 
flux by (|12p enables us to evaluate heat conductivity A according to Fourier's law. It's worth stressing that, as shown 
in the latter figures, unlike one-dimensional chains such as the one discussed above, fully three-dimensional models 
do predict normal heat conduction even when using harmonic potentials as ((5J), §§§ and (|10[) . Interestingly, in our 
simulations it is possible to drop out at will some of the interaction terms Vb, V a , V r b and V n b and investigate how 
temperature profile and thermal conductivity A are affected. It was found that potentials Vb and V a are strictly needed 
to avoid a collapse of the nanotube. Results corresponding to several setups are reported in Fig. 3 and Table 2. It is 
worth stressing that, for all simulations in a vacuum, nonbonded interactions V n b proved to have a negligible effect on 
both the slope of temperature profile and heat flux at steady state. On the contrary, the torsion potential V r b does 
have impact on the temperature profile while no significant effect on the heat flux was noticed: As a consequence, 
in the latter case, thermal conductivity shows significant dependence on V r b- More specifically, the higher torsion 
rigidity the flatter the temperature profile. 



IV. THERMAL BOUNDARY CONDUCTANCE OF A CARBON NANOFIN IN WATER 



A. Steady state simulations 



In this section, we investigate on the heat transfer between a carbon nanotube and a surrounding fluid (water) . The 
latter represents a first step towards a detailed study of a batch of single carbon nanotubes (or small bundles) utilized 
as carbon nanofins to enhance the heat transfer of a surface when transversally attached to it. To this end, and limited 
by the power of our current computational facilities, we consider a (5,5) SWNT (with a length L < 14 [nm]) placed 
in a box filled with water (typical setup is shown in Fig. 4). SWNT end temperatures are set at a fixed temperature 
That = 360 [A"], while the solvent is kept at T w = 300 [A"]. The carbon- water interaction is taken into account by means 
of a Lennard- Jones potential between the carbon and oxygen atoms with a parameterization (eco, °~Co) reported 
in table 1. Moreover, nonbonded interactions between the water molecules consist of both a Lennard- Jones term 
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FIG. 3: A (5,5) SWNT (green) is surrounded by water molecules (blue, red). Nose-Hoover thermostats with temperature 
That = 360 [-K] are coupled to the nanotube tips, while water is kept at a fixed temperature T w = 300 [-/£"]. After a sufficiently 
long time (here 15[ns]), a steady state condition is reached. MD simulation results (in terms of both temperature profile and 
heat flux) are consistent with a continuous one-dimensional model as described by eq. (|17[) and (|18[). Image obtained using 
VEGA ZZ 



between oxygen atoms (with too , &oo from table IIVB[) and a Coulomb potential: 

V c (r ij ) = J l-^, (13) 

where eo is the permittivity in a vacuum while qi and qj are the partial charges with qo — —0.82 e and qjj — 0.41 e 
(see also (H). 

We notice that, the latter is a classical problem of heat transfer (pictorially shown in Fig. 5), where a single 
fin (heated at the ends) is immersed in fluid maintained at a fixed temperature. This system can be conveniently 
treated using a continuous approach under the assumptions of homogeneous material, constant cross section S and 
one-dimensionality (no temperature gradients within a given cross section) (37} • In this case, both temperature field 
and heat flux only depend on the coordinate x, and the analytical solution of the energy conservation equation yields, 
at the steady state, the following relationship: 

f (x) = Me~ mx + Ne mx , (14) 

where T (x) = T (x) — T w denotes the difference between the local temperature at an arbitrary position x and the 
fixed fluid temperature T w . Let a and C be the thermal boundary conductance and the perimeter of the fin cross 
sections, respectively, m is linked to geometry and material properties as follows: 



whereas the two parameters M and TV are dictated by the boundary conditions, T (0) = T (L) = That (or cquivalently, 
due to symmetry, zero flux condition: dT/dx(L/2) = 0), namely: 

m.L/2 -mL/2 

M = f(0) T~F> FT?. N = f(0) — — . (16) 



Thus, the analytical solution (flUj) takes a more explicit form 

sh [i 

cosh {mL/2) 



~, , ~. . cosh \m (L/2 — x)] 

T(x)=T(0) \) ' /0 . h , (17) 



whereas the heat flux at one end of the fin reads: 



q = mXST (0) tanh (mL/2). 



(18) 
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\y 

FIG. 4: A single fin surrounded by a fluid can be studied by a one-dimensional continuous model, where all fields are assumed 
to vary only along the x-axis. 




FIG. 5: Steady state molecular dynamics (MD) simulations. Dimensionless temperature computed by MD (symbols) versus 
temperature profile predicted by continuous model (line), eq. (|17[) . Best fitting is achieved by choosing mL/2 — 0.28. Case 
with computational box 2.5 x 2.5 x 14[nm ]. 



In the setup illustrated in Fig. 4 and 5, periodic boundary conditions are applied in the x, y and z directions and all 
simulations are carried out with a fixed time step dt = l[/s] upon energy minimization. First of all, the whole system 
is led to thermal equilibrium at T = 300 by Nose-Hoover thermostatting implemented for 0.8 [ns] with a relaxation 
time tt = 0.1[ns]. Next, the simulation is continued for 15 [ns] where Nose-Hoover temperature coupling is applied 
only at the tips of the nanofm (here, the outermost 16 carbon atom rings at each end) with That = 360 [if], and 
water with T w = 300 [K] until, at the steady state, the temperature profile in Fig. 6 is developed. Moreover, pressure 
is set to l[6ar] by Parrincllo- Rahman prcssostat during both thermal equilibration and subsequent non-equilibrium 
computation. We notice that the above molecular dynamics results are in a good agreement with the continuous 
model for single fins if mL/2 = 0.28. Hence, this enables us to estimate the thermal boundary conductance a s t 
between SWNT and water with the help of eq. (p~5|) : 

m 2 XS 

a s t = (19) 

The thermal conductivity A has been independently computed by means of the technique illustrated in the sections 
above for the SWNT alone in a vacuum. Results for a nanofin with L = 14 [nm] are reported in Table 2. We stress that 
heat flux computed by time averaging of the Nose- Hoover parameter £ (see eq. (fl~2|)) is also in excellent agreement 
with the value predicted by the continuous model through eq. ([T5]) . For instance, with the above choice mL/2 = 0.28, 
for (5,5) SWNT with L = 10 [nm], L NH = 2[nm] in a box 5 x 5 x 14[nm 3 ] we have: - (f) N f k b T = 3.11 x 10~ 8 [VF] 
while 

q a = m\Sf (0) tanh (mL/2) = 3.14 x 10~ 8 . (20) 

We stress that L^h is the axial length of the outermost carbon atom rings coupled to a thermostat at each end of 
a nanotube. Finally, a useful parameter when studying fins is the thermal efficiency f2, expressing the ratio between 
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FIG. 6: Transient simulations. Temperature evolution as predicted by NVE molecular dynamics. Best fitting of exponential 
decay of the temperature difference Tcnt — T m is achieved by choosing: Td = 41 [ps]. 



the exchanged heat flux q and the ideal heat flux corresponding to an isothermal fin with T(x) = T(0), Vx G [0, L] 
[3^ ]. In our case, we find highly efficient nanofins: 

n = _± = m\Sf(Q) tanh(m£/2) = tanh(mL/2) = Q ^ 
qid a st CT(0)L/2 mL/2 



B. Transient simulations 



The value of thermal boundary conductance between water and a single wall carbon nanotube has been assessed by 
transient simulations as well. Results by the latter methodology are denoted as a tr in order to distinguish them from 
the same quantities (a st ) in the above section. Here, the nanotube was initially heated to a predetermined temperature 
Thot while water was kept at T w < T] lot (using in both cases Nose-Hoover thermostatting for 0.6[?is]). Next, an NVE 
molecular dynamics (ensemble where number of particle N, system volume V and energy E are conserved) were 
performed, where the entire system (SWNT plus water) was allowed to relax without any temperature and pressure 
coupling. Under the assumption of a uniform temperature field Tcnt{£) within the nanotube at any time instant t 
(Biot number Bi < 0.1), the above phenomenon can be modeled by an exponential decay of the temperature difference 
(Tcnt — T w ) in time, where the time constant Td depends on the nanotube heat capacity ct and the thermal heat 
conductance atr at the nanotube- water interface as follows: 

T d = — . (22) 
atr 

In our computations, following (39j . we considered the heat capacity per unit area of an atomic layer of graphite 
ct = 5.6 x 10 _4 [Jm _2 A' -1 ]. The values of and at r have been evaluated in different setups, and results arc reported 
in the table 2. It is worth stressing that values for thermal boundary conductance obtained in this study are consistent 
with both experimental and numerical results found by others for single wall carbon nanotubes within liquids [39l l4fj| . 



V. CONCLUSIONS 



In this work, we first investigate the thermal conductivity of single wall carbon nanotubes by means of classical 
non-equilibrium molecular dynamics using both simplified one-dimensional and fully three-dimensional models. Next, 
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TABLE I: Parameters for carbon-carbon, carbon- water and water- water interactions are chosen according to [13] and [34| . 



Carbon-Carbon interactions 


k b 


47890 kJmor^nm-' 1 


4, 


562.2 kJmor 1 




25.12 kJmor 1 


tec 


0.4396 kJmor 1 


VCC 


3.851 A 


Carbon-Oxygen interactions 


eco 


0.3126 kJmor 1 




3.19 A 


Oxygen-Oxygen interactions 


eoo 


0.6502 kJmor 1 


COO 


3.166 A 


Oxygen- hydrogen interactions 


<7o 


-0.82 e 


<7H 


0.41 e 



TABLE II: Summary of the results of MD simulations in this work. Single wall nanotubes with chirality (3,3), (5,5) and 
(15, 0) are considered, and several combination of interaction potentials are tested. In the first column, B, A, D and LJ stand 
for bond stretching, angular, dihedrals and Lennard- Jones potentials, respectively, while Bw denotes bond stretching with a 
smaller force constant k\j — 42000 \kJmol~ 1 K" 1 ] according to [25| . Simulations are carried out both in a vacuum (vac) and 
within water (sol). 



Chirality, Case 


Box 
[nm 3 ] 


Lnh 
[nm] 


L 
[nm] 


A 

1 m. K \ 


Otst 

La k\ 


Cttr 


Td 

[ps] 


mi/2 


(5,5), BAD-LJ (vac) 


12 x 12 x 12 


1.5 


5.5 


67 










(5,5), BwAD-LJ (vac) 


12 x 12 x 12 


1.5 


5.5 


64 










(5,5), BAD (vac) 


12 x 12 x 12 


1.5 


5.5 


65 










(5,5), BA (vac) 


12 x 12 x 12 


1.5 


5.5 


49 










(5,5), BwA (vac) 


12 x 12 x 12 


1.5 


5.5 


48.9 










(5,5), BAD-LJ (vac) 


20 x 20 x 20 


2 


10 


96.9 










(5,5), BAD-LJ (vac) 


105 x 105 x 105 


25 


25 


216.1 










(5,5), BAD-LJ (sol) 


2.5 x 2.5 x 14 


2 


10 




5.18 x 10' 






0.28 


(5,5), BAD-LJ (sol) 


4 x 4 x 14 


2 


10 




5.18 x 10' 






0.28 


(5,5), BAD-LJ (sol) 


4 x 4 x 14 





14 






1.70 x 10 ' 


33 




(5,5), BAD-LJ (sol) 


5x5x5 





3.7 






1.37 x 10 ' 


41 




(15,0), BAD-LJ (sol) 


5x5x5 





4.7 






1.60 x 10' 


35 




(15,0), BAD-LJ (sol) 


5x5x5 





3.8 






1.43 x 10' 


39 




(3,3), BAD-LJ (sol) 


5x5x5 





3.7 






8.90 x 10 b 


63 





based on the latter results, we have focused on the boundary conductance and thermal efficiency of single wall carbon 
nanotubes used as nanofins within water. More specifically, toward the end of computing the boundary conductance 
a, two different approaches have been implemented. First, a = a s t was estimated through a fitting procedure of 
results by steady state MD simulations and a simple one-dimensional continuous model. Second, cooling of SWNT 
(at Tqnt) within water (at T w ) was accomplished by NVE simulations. In the latter case, the time constant Td of the 
temperature difference (Tqnt — T w ) dynamics enables to compute a = at r - Numerical computations do predict pretty 
high thermal conductance at the interface (order of 10 7 [Wm^K -1 ]), which indeed makes carbon nanotubes ideal 
candidates for constructing nanofins. We should stress that, consistently with our results a s t > cttr, it is reasonable 
to expect that a s t represents the upper limit for the thermal boundary conductance, due to the fact that (in steady 
state simulations) water is forced by the thermostat to the lowest temperature at any time and any position in 
the computational box. Finally, it is useful to stress that, following the suggestion in [41|, all results of this work 
can be generalized to different fluids using standard nondimcnsionalization techniques, upon a substitution of the 
parameterization (ecd &Co) representing a different Lennard- Jones interaction between SWNT and fluid molecules. 
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VIII. METHODS 

The carbon nanotubes geometries simulated in this paper were generated using the program Tubegen [32j . while 
water molecules were introduced using the SPC/E model implemented by the genbox package available in GROMACS 
[3II ]. Numerical results in this work are based on non-equilibrium molecular dynamics where the all-atom forcchclds 
OPLS-AA is adopted for modeling atom interactions. Visualization of simulation trajectories is accomplished using 
VEGA ZZ H|. 
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